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Abstract. We study the rheology of semi-dilute solutions of the sodium salt of calf-thymus DNA in the 
linear and nonlinear regimes. The frequency response data can be fitted very well to the hybrid model 
with two dominant relaxation times r and n. The ratio ^ ~ 5 is seen to be fairly constant on changing 

the temperature from 20° C to 30° C. The shear rate dependence of viscosity can be fitted to the Carreau 
model. 



PACS. 87. 15. He Dynamics and conformational changes of biomolecules 
- 83.60Df Nonlinear viscoelasticity 



83.60Bc Linear viscoelasticity 



1 Introduction 

Deoxyribonucleic acid (DNA) is a key constituent of the 
nucleus of living cells and is composed of building blocks 
called nucleotides consisting of deoxyribose sugar, a phos- 
phate group and four nitrogenous bases - adenine, thymine, 
guanine and cytosine j[) . X-ray crystallography shows that 
a DNA molecule is shaped like a double helix, very much 
like a twisted ladder S. The ability of DNA to contain 
and transmit genetic information makes it a very impor- 
tant biopolymer that has been the subject of intense sci- 
entific research in recent years. DNA macromolecules are 
charged and depending on their molecular weight, the con- 
formation could vary between a rigid rod and a flexible 

" coil. The ability of flexible DNA molecules to twist, bend 
and change their conformation under tension or shear flow 
have been extensively studied both theoretically and ex- 

, perimentally J3|,|]]. The linear viscoelastic moduli of calf- 
thymus DNA solutions have been measured by Mason et 
al. at room temperature in the concentration range 1-10 

. mg/ml [|J. They have noted that the measured viscoelastic 
spectra do not fit the standard reptation model for flcxi- 

, ble polymers . Measurements of the nonlinear rheology 
of entangled T4 Bacteriophage DNA molecules @] show a 
plateau region in the stress a after an initial Newtonian 
regime at very low shear rates 7, similar to the flow curves 
seen in surfactant gels HP . 

The highly charged DNA macromolecule is a typical poly- 
electrolyte and the effects of the long-range and intra- 
chain Coulomb interactions on its structure and dynamics 
are significant. The structure of short DNA fragments in 
aqueous salt solutions have been studied using small-angle 
neutron scattering p0| . Direct mechanical measurements 
of the elasticity of single A-DNA molecules show devia- 
tions from the force curves predicted by the freely-jointed 
chain model jy], leading to the conclusion that DNA in 
solution has significant local curvature 1 12 . Measurements 



of relaxation times of a single DNA molecule manipulated 
using laser tweezers and observed by optical microscopy 
[[l3| show a qualitative agreement with the dynamic scal- 
ing predictions of the Zimm model [l^| . In recent years, 
experiments have been carried out on the electrohydro- 
dynamic instability observed in DNA solutions under the 
action of a strong electric field Q . Electric fields of the 
proper frequency and amplitude lead to the formation of 
islands of circulating molecules, with the islands arranging 
into a herring-bone formation. Recent simulation studies 
on short, supercoiled DNA chains show that DNA is a 
glassy system with numerous local energy minima under 
suitable conditions Q. 

In this paper, we discuss our recent results on the lin- 
ear and nonlinear rheology of semi-dilute solutions of the 
sodium salt of calf-thymus DNA. We fit our frequency re- 
sponse data to the hybrid model Q and the shear viscos- 
ity vs. shear rate data to the Carreau model |1S| ]. These 
models are essentially for dilute polymer solutions and 
the value of the radius of gyration R g of the DNA macro- 
molecules calculated from the fits to the hybrid model [[l7j 
are highly overestimated. We believe that the polydisper- 
sity of the DNA chains, the overlap of the macromolecules 
in the semi-dilute concentration regime and the electro- 
static interaction between the ionized groups also need to 
be considered for a more accurate theoretical description 
of our experimental results. 



2 The hybrid model theory 

The hybrid model theory incorporating Zimm dynamics 
has been used in the dilute polyelectrolyte and polymer 
literature to explain the behaviour of the reduced elas- 
tic modulus G' r (lo) —\\Ta c ^Q G' (uj) and the reduced vis- 
cous modulus Gft(u>) —\im c ^o[G" (uj) — ujr] s ] observed in 
aqueous solutions of sodium poly(styrene sulphonate) | fl9| 
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and dilute aqueous samples of polymers such as separan 
AP-30, xanthan gum, carboxymethyl cellulose and Polyox 
WSR301 (20). Here, rj s is the solvent viscosity and c is 
the polymer concentration. The hybrid model is charac- 
terised by a series of relaxation times spaced according to 
the Zimm theory jl4j , together with one additional longer 
relaxation time r . The Zimm model treats a single poly- 
mer chain in the framework of a bead-spring model, with 
N beads connected by N+l springs, in the presence of 
hydrodynamic interactions. The hybrid model is charac- 
terised by a series of relaxation times spaced as in the 
Zimm theory |l4| , in addition to a longest relaxation time 
t . In the limit of infinite dilution, the intrinsic moduli for 
this model may be written as ]l7| 



G' R {w) = G oW 2 t 2 [1 + 



N 



1 + W 2 T 2 ( _£.) 



[1] 



and 3^=1.66 according 



= ^r [Go(l + u 2 T 2 )~ 1 +m 2 ] + dE*! LJTlCl l, 

[2], 

where r p = 4^, p = 1, 3, 5 
to the Zimm model. We recall that in the Zimm model, 
only the time scales r p , which correspond to internal mo- 
tions such as flexure of the chain, contribute to the stress 
relaxation process. The experimental data for the elas- 
tic and viscous moduli of poly(2-vinyl pyridene) pl| have 
been fitted to the predictions of the Zimm model. The 
peaks in the viscoelastic moduli observed as a function 
of polyelectrolyte concentration have been explained in 
terms of electrostatically driven polymer coil expansion 
and contraction. Computer simulations have shown that 
on increasing the charge on the macromolecule, the longest 
relaxation time t\ in the Zimm spectrum 22 is enhanced. 
The hybrid model was proposed by Warren et al. to ex- 
plain the infinite-dilution viscoelastic properties of the he- 
lical molecule PBLG (poly-7- benzyl -L -glutamate) ||l7f . 
The value of r for infinitely dilute solutions of PBLG, 
obtained from the fits of G' R and G R to the Eqns. 1 and 2 
respectively, has been associated with the rotational dif- 
fusion of the macromolecule. The fitted value of T\ has 
been explained in terms of the flexural modes of vibration 
of the helix damped by solvent viscosity. Okamoto et al. 
have explained the viscoelastic properties of dilute solu- 
tions of the polyelectrolytes poly (acrylic acid) and poly 
(methacrylic acid) using this model, but without extrap- 
olation to infinite dilution |2^,Q. The authors conclude 
that the relaxation of these polyelectrolytes involves the 
rotation of the whole macromolecule, together with inter- 
nal configurational changes. 



3 Experiment 

We have carried out rheometric measurements on semi- 
dilute solutions of the sodium salt of calf-thymus DNA, 
dissolved in Tris-EDTA buffer (pH adjusted to physiolog- 
ical conditions of 7.9) at a concentration of lmg/ml (over- 
lap concentration of calf-thymus DNA is c*=0.35 mg/ml 
for a good solvent). The DNA was purchased from Sigma 



Chemicals, India in lyophilized form. Calf-thymus DNA 
is a linear, double-stranded macromolecule, and consists 
of ~ 13,000 base pairs. Since 1 nucleotide = 324.5 Da, 
the estimated molecular weight of calf-thymus DNA is ~ 
8.4 xlO 7 Da. We have dispersed the lyophilized DNA pow- 
der in an aqueous buffer consisting of lOmM Tris, lOOmM 
NaCl, 50mM each of NaCl and KC1 and 5 mM MgCl 2 . The 
samples prepared in this way were allowed to equilibrate 
for a day at 4°C to inhibit degradation. The rheological 
measurements have been conducted at 20°C, 25°C and 
30° C in an AR-1000N Rheolyst stress controlled rheome- 
ter (T. A. Instruments, U. K.) equipped with temperature 
control and software for shear rate control. We have used 
a stainless steel cone-and-plate assembly with a diameter 
of 4 cm and an angle of 1°59' as the shearing geometry. All 
the experiments reported below have been carried out on a 
sample of concentration 1 mg/ml. Because all experiments 
have been performed on the same sample, we have allowed 
sufficient time between runs to ensure that the sample re- 
laxes fully to its equilibrium state before the start of the 
subsequent experiment. For the oscillatory experiments, 
the rheometer has the provision for showing the waveform 
of the response on the application of an oscillatory stress. 
Care has been taken to ensure a distortion-free response 
to the applied oscillatory stresses for all the experiments. 



4 Results 

In this section, we discuss our experimental results on the 
linear and nonlinear rheology of calf-thymus DNA (con- 
centration c=lmg/ml) at 20°C, 25°C and 30°C. 



4.1 Linear rheology 

Figs. 1, 2 and 3 show the frequency response data (i.e. 
the measured G'(u>) and [G"(u>) — lot] s ] vs. the angular 
frequency u, where the solvent viscosity r/ s = lxlO -3 Pa) 
at 20°C, 25°C and 30°C and the corresponding fits (shown 
by solid lines) to the hybrid model with p =1 (Eqns. 
1 and 2) over almost three decades of angular fre- 
quency. Insets of figs. 1, 2 and 3 show the calculation of 
the oscillatory stress a osc lying in the linear regime. Be- 
cause G'(oj) and G"(cu) do not change appreciably over 
the range 8mPa to 30mPa at an oscillatory frequency li 
of 0.628 rads -1 , we have controlled the oscillatory stress 
for subsequent frequency response experiments at <r osc = 
25mPa. The values of the fitting parameters G , r , G\, 
T\ and TO2 are shown in Table 1. The measured G'(u>) and 
[G"(u) — urjs] could not be fitted to the Tanaka model [^5| 
for rigid rods given by G' R (ui) — G w 2 r 2 (l +w 2 t 2 ) _1 and 
G R (u) — ujt [G (1 + cj 2 t 2 ) -1 + m 2 ] over the entire fre- 
quency range for the frequency response data correspond- 
ing to T=20°C, 25°C and 30°C [§|. It has been pointed 
out earlier by Mason et al. Q that the frequency response 
data for calf-thymus DNA could not be fitted to the stan- 
dard models of reptation dynamics. We also find that the 
data for G'(ui) and G"(uS) (£6| cannot be fitted to the Doi- 
Edwards model P,E7| for flexible polymers. 
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Fig. 1. Elastic modulus G' (uj) (open circles) and viscous mod- 
ulus [G"(uj) — u>rj 3 ] (open triangles) vs. angular frequency uj at 
T=20°C and a osc = 0.025Pa. The solid lines show the fits to 
the hybrid model The inset shows the plot of G'(uj) and 
G" (uj) vs. the oscillatory stress a OS c at uj = 0.628 rads -1 . 




co (rads 1 ) 

Fig. 2. Elastic modulus G' (uj) (open circles) and the viscous 
modulus [G"(uj) — uir/s] (open triangles) vs. angular frequency 
uj at T=25°C and a osc = 0.025Pa. The solid lines show the fits 
to the hybrid model The inset shows the plot of G '(uj) 
and G" (ui) vs. the oscillatory stress a osc at uj — 0.628 rads -1 . 

In fig. 4, we have plotted the dynamic viscosity rj*(uj) 
vs. uj, where rf(u) = (G" 2 +G" 2 ) a5 at 2 o°C (up-triangles), 
25°C (circles) and 30°C (down-triangles). 
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Fig. 3. Elastic modulus G'(uj) (open circles) and viscous mod- 
ulus [G" (uj) — ujrjs] (open triangles) vs. angular frequency uj at 
T=30°C and a osc = 0.025Pa. The solid lines show the fits to 
the hybrid model The inset shows the plot of G '(uj) and 
G" (uj) vs. the oscillatory stress aw at uj = 0.628 rads -1 . 

Table 1 . Values of the fitting parameters G a , t c , Gi , ri and m,2 
obtained by fitting the hybrid model to G' (uj) and [G"(uj)~ uir/a] 
respectively at T=20°C, 25° C and 30° C. 



T °C 




To(s) 


Go (Pa) 


Ms) 


Gi(Pa) 


m2(Pa) 


20 


G' 


2.33 


0.021 


0.430 


0.080 




20 


g" 


1.75 


0.027 


0.330 


0.096 


0.326 


25 


G' 


1.15 


0.0344 


0.246 


0.097 




25 


g" 


1.140 


0.0386 


0.234 


0.1028 


0.383 


30 


G' 


1.31 


0.021 


0.230 


0.062 




30 


g" 


1.14 


0.020 


0.22 


0.069 


0.335 



Table 2. Values of the fitting parameters r) , tr and m ob- 
tained by fitting the Carreau model |l8| to 77(7) vs. 7 at 
T=20°C, 25° C and 30° C, respectively. 



T °C 


rj (Pa-s) 


tr (s) 


m 


20 


0.060±0.002 


0.499±0.145 


0.2425 


25 


0.053±0.002 


0.34±0.11 


0.2422 


30 


0.041±0.08 


0.36±0.08 


0.2430 



4.2 Nonlinear rheology 

Fig 5 shows the plot of stress a versus shear rate 7 at 20° C 
(up triangles), 25°C (circles) and 30°C (down triangles) 
respectively. Initially upto 7 ~ 2s -1 , the slope of the flow 
curve is very close to the Newtonian value of 1. At 7 > 
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4. Dynamic viscosity rf (ui) vs. angular frequency u 
=20° C (up-triangles) , 25° C (circles) and 30° C (down- 
;les) . 
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Fig. 5. Plot of the stress a vs. shear rate 7 at T=20°C (denoted 
by up-triangles), 25° C (circles) and 30° C (down triangles). The 
solid lines show the linear fits to the plots in the two regimes 
given by 7 < 2 and 7 >2. Note that at 7 < 2, the slope of the 
fit is very close to 1. 



2s , shear-thinning occurs and the fits to a ~ j a give a 
— 0.55. In fig. 6, we have plotted the shear viscosity 77(7) 
versus shear rate 7 at (a) 20°C, (b) 25°C and (c) 30°C. 
The plots of viscosity are found to fit well to the Car- 
reau model |D| written as 77(7) = ^ 1+ P 2 ° T -i ya ■ The values 

of the parameters obtained by fitting the shear viscosity 
data to the Carreau model have been listed in Table 2. 
Note that the values of 770 obtained from these fits are 
found to decrease with increasing temperature. However, 
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Fig. 6. Plot of the shear viscosity 77(7) vs. shear rate 7 at (a) 
T=20°C, (b) 25°C and (c) 30°C. The solid lines are the fits to 
the Carreau model. 



the temperature window is not wide enough to ascertain 
a possible Arrhenius dependence of the viscosities. 



5 Discussions of the results 

Fig. 7 shows the plots of the fitting parameters G , G\, r , 
77 and 7712 obtained by fitting the frequency response data 
at T=20°C, 25°C and 30°C to the hybrid model (Eqns. 1 
and 2) vs. the temperature. It is seen that the relaxation 
times t and t\ decrease monotonically on increasing tem- 
perature, while the percentage changes in G and G\ are 
less significant than the changes in t and t\ (Table 1). 
r may be associated with the rotational diffusion of each 
DNA macromolecule as a whole, whereas the time scale 
Ti , which corresponds to the internal degrees of freedom of 
DNA, might signify the time scale of the longest flexural 
mode of the macromolecule. The ratio of the two relax- 
ation times ^ ~ 5 is found to be fairly independent of 
the temperature. In an experimental study on the heli- 
cal macromolecule PBLG [17|, it has been shown that the 
relaxation spectrum, which may be fitted to the hybrid 
model, consists of a dominant relaxation time t describ- 
ing the end-over-end rotation of the macromolecule, and 
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Fig. 7. Plot of the fitting parameters Go 
obtained from the fits to the hybrid model | 
response data at T=20°C, 



Gi, To, Ti and m,2 
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25° C and 30° C. The ratio ^, which 
is also plotted, shows a negligible variation with temperature. 



a time t\ associated with its flexural bending mode. As- 
suming calf-thymus DNA to be a flexible coil undergo- 
ing predominantly rotational diffusion, its radius of gy- 



ration estimated from the relation r 



0.325 x 6 = 



k B T 

]27[ is found to be 1.6/xm at 25°C. This value overesti- 
mates the predicted value of R g ~ 0.4/xm for calf-thymus 
DNA @. For a rigid polymer, the rotational relaxation 

time r may be written as r = ^fc^T [ m (^)l H i where 
L is the length of the polymer and d is its diameter. Tak- 
ing 24=100, a rotational relaxation time of 1 s at 25°C 
gives an estimate of the length L ~ 2/im. Owing to the 
large configurational entropy of DNA and the fact that 
the DNA is in the semi-dilute concentration regime, the 
possibility of DNA macromolcculcs existing as rigid rods 
is extremely remote. This is further confirmed by the poor 
fit of the frequency response data to the Tanaka model for 



rigid rods 1 25 . The reason behind the discrepancy of the 
calculated R g with other estimates || can be due to the 
application of the hybrid model to a finite concentration 
of DNA, this model being ideally applicable only in cases 
of infinite dilution. Clearly, the effects of the overlap of the 
DNA macromolecules need to be incorporated in the the- 



ory. The Rouse model |27j rather than the Zimm model 
might be a more suitable model to describe the results. 
We recall that in the Rouse model, the time-scales in the 
relaxation spectrum may be expressed as t p ~ p-, where 
p =1, 3, 5 etc. The polydispersity of the DNA chains and 
the electrostatic interaction between the ionized surface 
groups also need to be considered for an accurate descrip- 
tion of our results. In spite of these shortcomings in our 
explanation of t and t\ , we note that the frequency re- 
sponse data fits remarkably well to the hybrid model over 
almost three decades of angular frequency. 
Dilute solutions of rod-like macromolecules are known to 
show a power-law shear-thinning after an initial Newto- 
nian regime at low 7 20 , 2{| . The viscosity data can be 
fitted to the Carreau's model jl8), with the fitted values 
of ?7o increasing with decrease in temperature. Power-law 
shear-thinning in nonlinear flow experiments and Zimm 
dynamics in the linear rheology data has been previously 
observed in dilute solutions of polymers such as Separan 
AP30, xanthan gum etc. by Tarn et al. p0|. 



6 Conclusions 

We have studied the linear and nonlinear rheology of semi- 
dilute solutions of calf-thymus DNA. The frequency re- 
sponse data show excellent fits to the hybrid model pro- 
posed by Warren et al. O] . In order to understand clearly 
the physical significance of r and n, we need to ex- 
tend the results of the hybrid model to the case of semi- 
dilute polylectrolytes. The electrostatic repulsion between 
ionisable groups and the polydispersity of DNA macro- 
molecules need to be accounted for. The flow curves at all 
temperatures may be divided into two regimes: an initial 
Newtonian regime at 7 < 2s _1 followed by a region of 
shear thinning at 7 > 2s _1 , as indicated by a — 0.55 in 
the fits to a ~ <-y - 55 at all temperatures. The shear vis- 
cosities 77(7) measured in these experiments and plotted 
vs. 7 can be fitted to the Carreau model pl|. It will be 
interesting to study the rheology of this system as a func- 
tion of DNA concentration, in order to understand the 
changes in the rheological parameters across the overlap 
concentration c* . 
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